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1. Introduction 



The conformational properties of linear polymers near an attractive wall are well 
understood by now |T]]. The fact that conformal invariance concepts are applicable 
to the problem has been extremely helpful, especially in two dimensions (in which 
case the "wall" is a line) where these tools provide a number of exact values for 
critical exponents. In contrast, for branched polymers it has been shown that the 
underlying field theory is not conformal [^. Exact results on bulk properties of 
randomly branched polymers have, however, been obtained through a connection with 
the theory of Yang-Lee edge singularities [Q. The corresponding extension towards 
surface properties has been accomplished only recently |^ yielding, among others, the 
interesting prediction that the crossover exponent at the adsorption point has the 
hyperuniversal (dimension-independent) value = 1/2. This result applies for an 
impenetrable wall; penetrable surfaces have not been considered 0. The case in favour 
of hyperuniversality is built from the following elements ||^: {€) an exact calculation 
in c? = 3, by means of a correspondence between branched polymers and an epidemic 
process plus a supersymmetric mapping of the latter onto a semi-infinite one-dimensional 
Yang-Lee edge problem; (u) conformal invariance properties of the two-dimensional 
Yang-Lee problem leading to information on four- dimensional branched polymers near 
a surface; {iii) mean-field theory, expected to be valid for (i > 8; and {iv) perturbation 
theory in c? = 8 — e dimensions, all four of which yield cf) = 1/2. 

In the present work we use finite-size scaling and phenomenological 
renormalisation ^ ideas to study surface properties of site lattice animals (which 
provide a model for randomly branched polymers in a good solvent) in two dimensions. 
To this end, the correlation length for animals on infinite strips with free boundary 
conditions (FBC) is numerically calculated by diagonalisation of the corresponding 
transfer matrix p]. By imposing FBC one is enabled to assess surface critical behaviour, 
in particular the adsorption transition. Our main goal is to check on the proposed 
hyperuniversal value = 1/2 at the adsorption point. Accordingly, only impenetrable 
surfaces are considered throughout this paper. We start by applying standard one- 
parameter phenomenological renormalisation (PR) 0, reobtaining bulk quantities such 
as the critical fugacity Xc and the temperature-like exponent y = l/u which are 
very accurately known 0, ^. This is important as a check of the overall reliability 
of our procedures. We then search for a surface-driven transition, by introducing a 
distinct fugacity for occupied sites along the strip boundaries. A two-parameter PR 
analysis is carried out, by comparing correlation lengths on three strips of consecutive 



widths [10, |TT], Two non-trivial fixed points are found, which are respectively related 
to the ordinary (bulk-dominated) and adsorption transitions. The corresponding finite- 
size estimates of critical parameters and exponents are extrapolated. Our main final 
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result is (f) = 0.505 ± 0.015, providing support to the hyperuniversality conjecture P 
= 1/2 at least in two dimensions. 

2. Model and calculational procedure 

The generating function of the model, defined on a semi-infinite lattice, can be written 
as: 

^ = E Cn,n.x''x^' , (1) 

N,Ns 

where C^^n^ is the number of different configurations that can be built with a total of 
sites constrained to form one cluster, of which Ns are at the surface; x is the fugacity 
for site occupation and Xg = exp tg/ksT where — is the extra energy assigned to each 



site at the surface. It is expected on general grounds [T^ that the critical fugacity Xc 
will be a constant as a function of Xg, from small Xg up to the adsorption threshold 
given by some > 1. Upon approach to a point {x^,x^g) on the critical line in {x,Xs) 
space, the bulk correlation length ^ diverges as ^ ~ 6~'^ where the scaling field 5 is a 
suitable linear combination of Sx = x — x^ and 6xs = Xg — x^ . Close to the adsorption 
point, a second length ^g diverges with a different exponent Ug = 1/yg (and a different 
combination 6g of the variables 6x and 6xg). Physically, ^g measures the thickness of the 
adsorbed polymer layer. Thus the average number of surface contacts < Ng > scales 
asymptotically with the average number of sites < > as 



<Ng>r^<N>\ = yg/y . (2) 

Below the adsorption threshold one has = 0, while in the adsorbed phase = 1. At 
the threshold is expected to take on a non-trivial value. 

The correlation length ^L^x^Xg) along a strip depends on the largest eigenvalue 
A° {x,Xg) of the transfer matrix via : ^^(a;, x^) = — l/lnA^, {x,Xg). For large L and 
close enough to criticality of the corresponding (semi-) infinite system, the correlation 
length must scale, in terms of 5 and 5g defined above, as 

L-^U^^5g) = F{Ly5,Ly^5g). (3) 



Upon rescaling of one expects two non-trivial fixed points [|TT], |T^, |T3|: {i) the 
ordinary fixed point, which governs the critical behaviour of the unbound (bulk) phase, 
at which the surface interactions are irrelevant and thus exhibits yg < 0; and {ii) the 
special fixed point, corresponding to the adsorption transition, with yg > 0. 

We use strips of width L < 10 sites, both for square and triangular lattices. Building 
up the transfer matrix involves the analysis of connectivity properties of adjacent 
columns of occupied and empty sites. For the present case of animals on strips with 
FBC, this is a straightforward extension of earlier work on percolation and animals with 
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periodic boundary conditions (PBC) |^ and percolation with FBC JT^. The resulting 
matrix is rather sparse, owing to restrictions imposed by connectivity [Q: for L = 10 
on the square lattice, for instance, only 4.2% of the possible combinations of adjacent 
column states are allowed. 

Extrapolation of finite-width results must be dealt with carefully, especially as 
convergence of estimates produced with FBC is usually slower than that of their 



counterparts generated with use of PBC |T^, [T^, |T^ . In the present work, extrapolations 
toward L ^ oo have been done using the Bulirsch-Stoer (BST) algorithm []T7|, |T^. As 



extensively discussed elsewhere whenever the leading correction-to-scaling exponent 
uj is not known a priori BST extrapolations rely on keeping it as a free parameter within 
an interval guessed to be reasonable. Central estimates and error bars are evaluated 
self-consistently by selecting the range of u for which overall fluctuations are minimised. 

3. Results 

3.1. One-parameter renormalization 

We first consider no surface binding {xg = 1). We can then implement standard, one- 
parameter, PR in the usual way by looking for a finite-size estimate of Xc given by the 
fixed point a;^ of the implicit recursion relation: 

L L- 1 ' 

At the fixed point, an approximation to the bulk exponent y = 1/z/ is evaluated by ||^: 

\n{{d^L/dx)^*/{d^L-i/dx)a,i} 
= ML/L-l) 

In order to check on universality of critical amplitudes ||19|, we also calculate the 



(4) 



quantity Al = 2L/tx^l{x*^). Note that for a triangular lattice with FBC the strip width 
L = where is the number of sites across the strip. If the underlying field 

theory were conformal at the critical point, this would be an estimate of the exponent 
describing the decay of critical correlations along the surface, rjg. 

Our results are shown in table 1. Overall agreement with expected values, where 
these are available, is rather good. Universality of critical correlation-length amplitudes 
is satisfied within error bars. However, finite-size data show that the amplitude of 
corrections is much larger than for the corresponding cases of PBC (see e.g. table I 
of reference 0). Partially as a consequence of this, our extrapolated estimates for Xc 
and y are somewhat less accurate than those obtained with PBC 0. A second source 
of imprecision emerges when one considers the broad ranges allowed for the correction- 
to-scaling exponent u in table 1. Though some degree of subjectivity is inevitable 
when dealing with error estimation within the BST scheme, our results for u reflect 
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Table 1. Results from one-parameter PR. Uncertainties in last quoted digits are shown 
in parentheses. Extrapolations obtained by BST algorithm with correction-to-scaling 
exponent w in ranges shown. Expected values from Ref. 



Square Triangular 



L x*j^ yL Al xl UL Al 

3 0.298906 1.43406 1.02271 0.247186 1.39519 0.998368 

4 0.274596 1.45806 1.27009 0.221892 1.42363 1.25590 

5 0.263725 1.47490 1.44619 0.210650 1.44332 1.43950 

6 0.257947 1.48691 1.57871 0.204743 1.45791 1.57633 

7 0.254533 1.49588 1.68219 0.201287 1.46922 1.68213 

8 0.252364 1.50284 1.76528 0.199107 1.47828 1.76637 

9 0.250909 1.50841 1.83350 0.197652 1.48570 1.83506 

10 0.249890 1.51296 1.89052 0.196638 1.49190 1.89217 
Expected 0.246150(10) 1.5607(4) — 0.192925(10) 1.5607(4) — 
Extrapolated 0.2460(2) 1.55(1) 2.45(4) 0.1928(2) 1.55(1) 2.4(1) 
w 1.5(4) 1.5(5) 1.45(20) 1.5(4) 1.5(5) 1.5(5) 



the fact that, roughly for u between 1 and 2 the fluctuation estimates for fixed u 
(based on the spread between next-to- highest order estimates [Q) keep to the same 



order of magnitude. On the other hand, outside this interval fluctuations increase, and 
estimates deteriorate, quickly. This is to be compared e.g. to similar extrapolations 
for percolation with FBC [jl4| where usually one can pinpoint a much narrower band of 
values of uj within which fluctuations are minimised. 

3.2. Two-parameter renormalisation 

We next allow the surface interaction to vary. Similarly e.g. to studies of linear polymer 



adsorption [|Tl|, |T2|, an extra energy — is introduced for sites on either strip boundary, 
so that Xs = expes/ksT. L-dependent fixed points are obtained by comparing 

correlation lengths on three strips ||T0[| : 

^l{x*,x*) ^ ^L-i{x*,x*) ^ ^L-2{X*,X*) 

L L-1 L-2 ■ ^ ' 

In the present case these equations give two fixed points: the ordinary fixed point, 
which describes the behaviour of the unbound animal and the special fixed point which 
describes the animal's behaviour at the binding transition. Linearizing around the 
fixed points, the exponents y and i/s can be found from suitable partial derivatives 
evaluated at the fixed point in question []TU|, |r^, [T2[. Again we calculate the quantity 
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Al = '^L/'kE,l{x* 1 X*). Our results for the ordinary and special fixed points are displayed 
respectively in tables 2 and 3. 

As a general rule, finite-size estimates differ from their limiting {L oo) values 
by much smaller amounts than was the case in one-parameter PR. In several instances, 
though, convergence turns out not to be monotonic. Further, within the BST scheme 
we frequently find the following as the trial value of uo is increased from 0.45 to, say, 
6: (i) fiuctuation estimates at fixed u always decrease, and {ii) last-order approximants 
Qioj) vary monotonically, and seem to be converging towards fixed points(that is, 
dQ{uj)/duj —>■ 0). This is consistent with what is found from three-point extrapolations 
adjusting u for the best straight- line fit of data against L^'^ [^]: as a rule, u tends to 
converge to unrealistically high values. Thus, although strictly speaking there are no 



regions where the BST algorithm is stable with respect to to [0 in such cases, one can 
produce reasonably reliable estimates by looking at trends followed upon increasing uj. 
For the entries in tables 2 and 3 to which this picture applies we display the ranges of 
variation of last-order approximants corresponding to uj > uJmini with uJmin as given in 
the respective entry. 



3.2.1. The ordinary transition. For the ordinary transition, the corresponding fixed 
point can be found only for L > 6. The exact result i/^ = — 1 is expected to hold, as it is 



based on general properties of the ordinary transition of two-dimensional systems |20 
For both square and triangular lattices, extrapolations were performed discarding data 
for L = 6. Though, for the latter, these do not usually imply non-monotonic variation 
along the sequence, their inclusion would increase the scatter of extrapolates by at least 
one order of magnitude. For the non-universal on the triangular lattice, we have found 
neither an optimal range for tu, nor the smooth decrease of error as uj increases, described 
above. Thus we quote for an average of last-order estimates for 1.0 < oj < 4.0. In 
general, the final results for the ordinary fixed point show agreement to within less than 
0.5% with those of Ref. for the exact result Us = —1 fluctuations are higher, but 
still kept smaller than 2%. Universahty of the amplitude A is satisfied within 0.05%. 



3.2.2. The special transition. For the special fixed point on the square lattice we have 
discarded L = 5 and 6 data for x*, x* and Al on account of non- uniform convergence; 
otherwise, all data in table 3 have been used in extrapolations. For t/s on the square 
lattice, fiuctuations were approximately constant and small throughout the range of uj 
explored, so we quote an average of last-order estimates for 1.0 < uj < 4.0. While 
extrapolates from square-lattice results undoubtedly suffer as a result of the above- 
mentioned difficulties, application of the BST algorithm nevertheless gives a fairly 
accurate numerical picture. 

On the other hand, results for the triangular lattice fall into smooth, well-behaved 
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Table 2. Results from two-parameter PR at the ordinary fixed point. Uncertainties in 
last quoted digits are shown in parentheses. Extrapolations obtained by BST algorithm 
with correction-to-scaling exponent lu in ranges shown. 

(a) Square 



T 

Li 


rr-* 

X 




y 


Vs 


A 


6 


0.246282 


0.223963 


1.57916 


- 1.02151 


2.49768 


7 


0.246420 


0.232356 


1.57746 


- 1.05163 


2.48955 


8 


0.246352 


0.225924 


1.57451 


- 1.04907 


2.49477 


9 


0.246294 


0.217824 


1.57223 


- 1.04172 


2.50040 


10 


0.246254 


0.209961 


1.57057 


- 1.03480 


2.50516 


Expected 


0.246150(10)'" 





1.5607(4)=" 


- l'' 





Extrapolated 


0.24617(3) 


0.18(1) 


1.566(1) 


- 1.014(6) 


2.519(4) 




> 3 


■> 3 n 


> 3 


> 3 n 


> 3 






(b) Trianj 


^ular 






L 


X* 




y 


Vs 


Al 


6 


0.194510 


0.352877 


1.53424 


- 0.836625 


2.381520 


7 


0.193555 


0.289360 


1.54238 


- 0.804216 


2.445904 


8 


0.193238 


0.257802 


1.54615 


- 0.841891 


2.472772 


9 


0.193115 


0.240468 


1.54833 


- 0.880553 


2.485607 


10 


0.193056 


0.229288 


1.54982 


- 0.910094 


2.492939 


Expected 


0.192925(10)^ 




1.5607(4)=" 


- l'^ 




Extrapolated 


0.19296(1) 


0.185(15) 


1.555(2) 


~ 1.00(1) 


2.512(2) 




> 3.0 


c 


2.4(9) 


3.3(5) 


2.4(6) 



Ref. @ 
■Ref. 13 



■^No optimal uo found (see text). 

sequences from which we have extracted a set of very precise extrapolates. Our central 
estimate for y is higher by 0.4% than that of Ref. [Q, with non-overlapping error bars. 
Recalling that our error bars reflect uncertainties in the extrapolation procedure itself, 
and do not take into account systematic errors in the original sequence of finite-size 
results, we do not take this fact as necessarily meaning that our estimates conflict. 
Indeed, other instances are known |12| in which extrapolates from two-parameter 
PR differ slightly e.g. from exact results. 

4. Discussion and conclusions 



It can be seen from tables 2 and 3 that universality of critical amplitudes is satisfied 
to within error bars. The qualitative behaviour as one spans the distinct possibilities 
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Table 3. Results from two-parameter PR at the special fixed point. Uncertainties in 
last quoted digits are shown in parentheses. Extrapolations obtained by BST algorithm 
with correction-to-scaling exponent lu in ranges shown. 

(a) Square 



L 


X* 




y 


Vs 


Al 


5 


0.244202 


2.34075 


1.54294 


0.717355 


- 0.0980272 


6 


0.246045 


2.28235 


1.55829 


0.796012 


- 0.0579771 


7 


0.246033 


2.28278 


1.56222 


0.796973 


- 0.0583358 


8 


0.246088 


2.28049 


1.56514 


0.800373 


- 0.0561538 


9 


0.246108 


2.27951 


1.56697 


0.801618 


- 0.0551068 


10 


0.246123 


2.27877 


1.56825 


0.802352 


- 0.0542317 


Expected 


0.246150(10)'^ 


— 


1.5607(4)^^ 


— 


— 


Extrapolated 


0.24615(1) 


2.277(1) 


1.571(2) 


0.804(1) 


- 0.054(2) 


LU 


> 3.0 


> 3.0 


> 3.9 


b 


> 3.0 






(b) Trian 


gular 






L 


X* 




y 


Vs 




5 


0.190915 


2.85122 


1.53472 


0.701013 


- 0.117316 


6 


0.192174 


2.78810 


1.54574 


0.745358 


- 0.0887032 


7 


0.192573 


2.76501 


1.55200 


0.764248 


- 0.0761753 


8 


0.192735 


2.75442 


1.55608 


0.773683 


- 0.0695443 


9 


0.192813 


2.74873 


1.55889 


0.778886 


~ 0.0655337 


10 


0.192856 


2.74527 


1.56086 


0.781982 


- 0.0628311 


Expected 


0.192925(10)^^ 




1.5607(4)^^ 






Extrapolated 


0.19292(1) 


2.745(5) 


1.5663(5) 


0.7888(5) 


- 0.054(3) 


LU 


3.5(6) 


2.7(4) 


3.6(6) 


3.5(7) 


2.0(4) 



-Ref. @ 

'^No optimal u; found (see text). 



is similar to that of the corresponding t] of hnear polymers. For PBC A — 0.68 p], 
with FBC, at the ordinary transition A ^ 2.51; at the special transition A ^ —0.054. 
For linear polymers rjbuik = 5/24; r/f'"'^ = 5/4; r^f^ = —1/12 |Tl|, |l^. Unfortunately the 
analogy does not seem to go beyond this level. 

The result = 2.277(1) for the adsorption threshold on the square lattice compares 
well with, and is more precise than, the series estimate 2.25(5) It would be 

interesting to check whether using our value of Xg in the series analysis would improve 
other results. 

Turning now to the crossover exponent = ys/y, the safest course seems to be 
separately extrapolating the sequences for yg and y, and then calculating the ratio of 
final estimates. From the square-lattice data of table 2 one would get = 0.511 ±0.002, 
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while for the triangular lattice of table 3 one gets (f) = 0.5035 ± 0.0005. Though, as 
mentioned earlier, the above error bars do not take into account systematic errors, it is 
desirable to have an estimate of these. In order to do so, we refer to the similar case of 
adsorption of two-dimensional linear polymers, where conformal invariance asserts that 
y = 4/3 and i/s = 2/3, thus 0=1/2 exactly; two-parameter PR sequences for a locally 
directed (but globally isotropic) square lattice extrapolate respectively to ~ 1.339 and 
ys ~ 0.679 |T2[, which gives ~ 0.507, just under 2% off the exact value. For fully 
isotropic lattices, the extrapolation for z/g overshoots the exact value by 4% |TT[. In 



this latter case the corresponding extrapolations of y have not been published; however, 
finite-lattice data point towards a value somewhat larger than 4/3 |Tl[], thus one expects 
at least a partial compensation of errors for the ratio ys/y- Assuming systematic errors of 
order 2% in the extrapolated exponents also for the present case, one gets = 0.51 ±0.01 
(square lattice) and 0.50 ± 0.01 (triangular). Systematic deviations are thus estimated 
to increase uncertainties by at least one order of magnitude over those coming from 
extrapolation procedures. Our final result must encompass both the latter central 
estimates and allow for the spread between them, plus their own inherent uncertainties. 
Thus we have = 0.505±0.015. This is consistent with the hyperuniversality conjecture 
0=1/2 1^, and to be compared with the series result = 0.6 ± 0.1 |^ . 



We have studied surface properties of randomly branched polymers in two 
dimensions. Our estimates of bulk quantities such as critical fugacity Xc and critical 
exponent y are in very good agreement with results obtained with PBC |^. We have 
checked that universality of critical amplitudes |TP|] holds in all instances investigated. 
The adsorption threshold for the square lattice has been located with greater accuracy 
than previously available . The crossover exponent at the adsorption point satisfies, 
within error bars, the recent hyperuniversality conjecture = 1/2 f^. 



Acknowledgments 

The author thanks M Henkel for interesting conversations, and Departamento de Fi'sica, 
PUC/RJ for use of their computational facilities. This research is supported by CNPq, 
FINEP and CAPES. 



References 



[1] Eiscnriegler E 1993 Polymers Near Surfaces (Singapore: World Scientific) 

[2] Cardy J L 1987 Phase Transitions and Critical Phenomena Vol 11 ed C Domb and J L Lebowitz 

(London: Academic) 
[3] Miller J D and De'Bell K 1993 J. Physique I 3 1717 
[4] Parisi G and Sourlas N 1981 Phys. Rev. Lett. 46 871 
[5] Janssen H K and Lyssy A 1994 Phys. Rev. E 50 3784 



10 

[6] Barber M N 1983 Phase Transitions and Critical Phenomena Vol 8 ed C Domb and J L Lebowitz 
(London: Academic) 

[7] Nightingale M P 1990 Finite Size Scaling and Numerical Simulations of Statistical Systems ed V 
Privman (Singapore: World Scientific) 

[8] Derrida B and de Seze L 1982 J. Physique 43 475 

[9] Derrida B and Stauffer D 1985 J. Physique 46 1623 
[10] Derrida B and Herrmann H J 1983 J. Physique 44 1365 
[11] Veal A R, Yeomans J M and Jug G 1991 J. Phys. A: Math. Gen. 24 827 
[12] de Queiroz S L A and Yeomans J M 1991 J. Phys. A: Math. Gen. 24 1874 

[13] Binder K 1983 Phase Transitions and Critical Phenomena Vol 8 ed C Domb and J L Lebowitz 

(London: Academic) 
[14] de Queiroz SLA 1995 J. Phys. A: Math. Gen. 28 L363 
[15] Derrida B and Vannimenus J 1980 J. Physique Lett 41 L473 
[16] Burkhardt T W and Guim I 1985 J. Phys. A: Math. Gen. 18 L25 
[17] Bulirsch R and Stoer J 1964 Numer. Math. 6 413 
[18] Henkel M and Schiitz G 1988 J. Phys. A: Math. Gen. 21 2617 
[19] Privman V and Fisher M E 1984 Phys. Rev. B 30 322 
[20] Burkhardt T W and Cardy J L 1987 J. Phys. A: Math. Gen. 20 L233 
[21] Foster D P and Seno F 1993 J. Phys. A: Math. Gen. 26 1299 



